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ANALYTIC ESTIMATE OF THE ORDER PARAMETER FOR MAGNETIC 
CHARGE CONDENSATION IN QCD 
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The order parameter for monopole condensation is computed in terms of gauge invariant field 
strength correlators. Important properties emerge of the correlators in the confined phase, which 
could not be extracted by existing numerical determinations on the lattice. 
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1. INTRODUCTION 

A possible mechanism for confinement of 
color is Dual Superconductivity of the Vac- 
uum 1 2 3 . The idea is that magnetic charges 
condense in the vacuum : the chromoelectric 
field acting between a q — q pair is channeled 
by dual Meissner effect into flux tubes with 
energy proportional to the distance. 

Flux tubes are indeed observed in lattice 
simulations , and the corresponding field con- 
figurations are consistent with expectations 

4 

An order parameter (/x) has been defined 
5 6 7 8 9 which j s vacuum expectation 
value of a magnetically charged operator \i. 

Numerical simulations show that (fi) =/= 
in the confined phase , signaling Higgs break- 
ing of magnetic U(l) symmetry , (fi) — in 
the deconfincd phase. Deconfincment is an 
order-disorder phase transition. 

For SU (N) gauge group there exist (A^ — 
1) magnetic charges 10 and the corresponding 
operators fi a , (a = 1, .., N — f ) can be written 
as 



fi a (x, t) = e 2 <> 



d 3 yb ± (x-y)Tr(^E)(y,t) 



(1) 



is the field of a Dirac 
= and V A b i = 4+ 



q is an integer b± 
monopole with V&x 
Dirac string. 

<S> a (x) = U(x,y)$ a d W(x,y) 

with U(x,y) an arbitrary gauge trans- 
formation which we shall take as a parallel 



transport to x from a reference point y along 
a path C , and 



< 



N -a N -a 



> . < -(N-a) > 



iV' N >[> 



n^diag{ N ,, N 

(/i a ) is gauge invariant. In the gauge in which 
$ a = $2 (Abelian Projection) 

H a (x,t) = eH/^M^itf.*) (3) 

where E\ is the transverse chromoelectric 
field along the color direction T a 

a . a + 1 

T a = diag(0 1 ...0,1 , -1,0, ..0) (4) 

and , since E a ± is the conjugate momentum 
to A a ± one has 

H a (5!,t)\A a ± (j?,t)) = \A±(#,t) + ±b x (x-y)) 

(5) 

(/Lt a ) is the ratio of two partition functions 7 
with the action proportional to /3 = ^ . At 
fi = \i a = 1. It proves convenient to use 
instead of (fi a ) the susceptibility 

„ d 



P 



83 



ln( f i a ) 



in terms of which 



(6) 



(7) 



p a is easily measured by lattice simulation, 
and this has been done for pure gauge theo- 
ries compact U(l) 7 , 5*7(2) 8 , 51/(3) 9 , and 
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for Nf — 2 QCD 11 . In all cases the result is 
the following, in the thcrmodynamical limit 
V = Ll -> oo 

1) T < T c (Confined ) 
p a — > finite as V — > oo , 
which by Eq(7) implies (^i a ) 7^ 

2) T > T c (Deconfined phase) 
p a w -|c|L s + c' or (/i a ) = 

3) T w T c (Critical region) 

■4«/(tL;) 
(finite size scaling) 

Here r = (1 — ^-) , z/ is the critical index 
of the correlation length of the order param- 
eter. 

p a is independent of the choice of the 
abclian projection 12 13 14 . 

From the definition of (p a ) it follows 



becomes 



■ f d 3 Vl f d 3 j/2*? lii2 {yi-V2)b'l (j/i)bi 2 (fc) 



(9) 



or, since (5 = -p- Eq(6) becomes 



P = — 



q 2 8 



yi 



J/2 



16N dp 1 

K,M-m)bl(yi)bl(m)} (">) 

We shall identify i2 with the two point 



3 ( — ) 



correlators defined with a straight line paral- 
lel transport which are measured on the lat- 
tice 171819 , and are used as an input in the 
stochastic model of QCD. In fact for those 
correlators (E a E b ) = 5 ab <£ , so that p a is in- 
dependent on a , and this agrees with lattice 
determinations of p a 9 . 

The validity of the cluster expansion is 
usually justified at large distances, and we 
are looking for infrared properties in the 
study of confinement. 
yn){($ a ■E)i 1 y 1 ($ a .E)i n yjfi) A direct check can of be obtained by 

studying the dependence on the monopolc 
charge q. The truncated p is proportional 
to q 2 : higher correlators would introduce 
terms proportional to higher powers of q. 
Old data 8 15 agree with proportionality to 



d z yi ..d z y n b l l{x - yi). 



notation 



We have used the simplified 
(<t> a E) hS = Trl<t> a (t,x)E*(t,x)}. 

The vev's in eq(8) are gauge invariant 
field strength correlators. We shall identify 
these correlators with those of the stochastic 
vacuum approach to QCD 20 21 22 . 



2. Cluster Expansion of (fi a ) 

The idea of the stochastic vacuum approach 
is to approximate systematically the mul- 
tiple gauge invariant field strength correla- 
tors by products of lowest order correlators 
in a systematic cluster expansion. The ex- 
pansion is then truncated and only clusters 
up to order two are retained. Since the 
one point cluster {(& a .E)) — 0, only two 
point clusters will survive the truncation, 
and hence only even terms in the expansion 
Eq(8) , which will be approximated as prod- 
ucts of two point correlators f (yi — y 2 ) = 
({<5> a .E) lu $ l (<5> a .E) l2 ^ 2 ). The combinatorial 
factor of the term 2n is (2n — 1)!! and Eq(8) 



q but a more systematic study of this issue 
is planned. 

3. The Field Correlators 

There exists a general parametrization of 
field strength correlators dictated by invari- 
ance arguments 20 21 

V{z 1 ,z 2 )F^. V2 {z 2 )V\z 1 ,z 2 )) (11) 

(D(z! - z 2 )[S lllll2 S VlV2 - <W 2 <W 2 ] 

1 d 

+ 7:^ — [Di(zi - z 2 )(z li2 8 UlU2 - z V2 5 Vlli2 ]+ 

2 Cz Ml 



1 d 



2dz p 



■[D(zi - z 2 )(^ 2 <W 2 ~ 



(12) 
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At T ^ the electric field correlators are not 
equal to those of the magnetic field, so one 
has four form factors, De, Die, Dh, Dm ■ 
For the electric field correlators [i\ = [i2 = 
0, v\ = ii, i>2 = 12 and 

^ b ll2 =5 ab [8 lll2 (DE + \D 1E ) + ^-] (13) 

In the convolution with b± the derivative 
terms give 0. For the same reason the re- 
placement can be performed 

X. . . X _ fc 'l fc '2 

Vl!l2 ' u lll2 fe2 

Going to the Fourier transform Eq(10) 
becomes 



P =- 



q 2 d 
16 9/3 



IP 



d 3 k 
(2tt) s 



D E (k 2 )j^(k 2 S knkz2 -k n ki 2 )} (14) 

where Z)g(fc 2 ) is the Fourier transform of 
(D E + \Di E ). The identity 

(k 2 d lJ -k n h 2 )b l l(k)b i l(-k) = \H(k)\ 2 (15) 

, together with the explicit form of H(k) 

H(k) = kAb ± (k) (16) 

with n the direction of the Dirac string ( we 
shall call it z), gives 



\H(k)\ 2 = -^ + ^ (17) 



and for p a 



a I 2 d ro f d3k I 1 1 \ 2m 

P = wt^J j^^-ki^ 

(18) 

Here f(k 2 ) = -^D E (k 2 ). Again we notice 
that identifying our correlators with those of 
the stochastic model implies that p a is inde- 
pendent both on a and on the abelian pro- 
jection. 

At large /3(deconfincd phase) f(k 2 ) can 
be approximated by first order perturbation 
theory 

1 



the only dependence on (3 is the explicit fac- 
tor in Eq(18) so that 

d 3 k 1 1 1 



r) (20) 



167V J (2ir) 3 k y k 2 k 2 ' 

The integral is easily evaluated with UV cut- 
off - (a the lattice spacing) and IR cut-off 
(L s the spacial lattice size). The result 

is 

^ = lfjV (L) 2 [-^)Ls + 2ln{L s )+con S t] 

^ (21) 

which, as explained in Sect. 1 means (^ a ) = 
in the thermodynamical limit L s — > oo 

The two point correlators have been 
measured on the lattice both at T — and 
at finite T. They are well represented in the 
range of distances .1/m < x < lfm by a 
form 17 



D, 



Ae x b + 



(22) 



with A a w 2Ab and Ab w .3/m. A and b 
are independent on (3 from T = up to 
T « .95T C Approaching further T c A de- 
creases rapidly to zero. 19 . 

In fact the above parametrization can- 
not be valid at shorter distances, where the 
operator product expansion and the non- 
existence of condensates of dimension less 
than 4 require that 

' : -]+c+dx 2 (23) 



2Nk 



(19) 



2 l {x + ie) 4 (x-ie) 4 ' 

The prescription on the singularity is dic- 
tated by the match to perturbation theory. 
At larger distances a stronger infrared cut- 
off at some distance A is expected since col- 
ored particles cannot propagate at infinite 
distance. This feature needs a further nu- 
merical investigation on the lattice. Up to 
T w .95T C the only dependence on in 
Eq(14) is again the explicit factor so that the 
result is the same as Eq(21) with the lattice 
size L s replaced by the infrared cut-off A 

^ = 16" (2^F ^ V^A + 2ln(K) + const] 

(24) 
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This expression gives a finite value of p a for 
T < T c and hence (p a ) ^ which means 
dual superconductivity for any finite value of 
the UV cut-off a. However in the continuum 
limit o^0p° diverges so that (/i a ) needs a 
renormalization. This is similar to what hap- 
pens for the Polyakov line in the quenched 
theory. 23 . Existing Lattice data support this 
statement [See Fig(2) of ref. 8 ] ,but a more 
systematic investigation is planned. When 
the critical temperature is approached both 
the IR cut-off A and the coefficient A in 
Eq(22) strongly depend on (3 :A diverges and 
A tends to zero. A more detailed calculation, 
which will be reported elsewhere gives 



g 2 d 
16N 9/3^ l (2tt) 



A 



+ const) + 2N XtA E (1 + — )] (25) 

ir\ b 

Numerical determinations of the tempera- 
ture dependence of Ae around T c exist. Not 
much is known about the behavior of A . Fur- 
ther study is needed to understand how the 
scaling law of p a described in Sect 1 works 
in Eq(25), and this is already on the way. 
In conclusion an interesting interrelation has 
been found between the large distance be- 
havior of correlators an confinement, which 
deserves further study. 
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